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Abstract The failure criterion is an essential part of all 
strength calculations of design. It was shown in the past that 
the tensor-polynomial equation could be regarded as a 
polynomial expansion of the real failure surface. Now it is 
shown that the third-degree polynomial is identical to the 
real failure criterion. It is also shown that the second-degree 
part of the polynomial is identical to the orthotropic exten¬ 
sion of the von Mises criterion for initial yield. The third- 
degree polynomial hardening terms of the criterion are also 
shown to incorporate the earlier theoretical explained 
mixed-mode I-II fracture equation, showing hardening to 
be based on hindered microcrack extension. For uniaxial 
loading, the failure criterion can be resolved in factors, 
leading to the derivation of extended Hankinson equations. 
This allows the relations between the constants of the total 
failure criterion to be elucidated, which is necessary for data 
fitting of this criterion and providing a simple method to 
determine the constants by the simple uniaxial, oblique- 
grain compression and tension tests. Based on this, the 
numerical failure criterion is given with the simple lower 
bound criterion for practice and for the codes. 

Key words Timber ■ Failure criterion • Tensor-polynomial • 
Extended orthogonal von Mises criterion ■ Extended 
Hankinson equations 


Introduction 

Design and control calculations are normally based on limit 
analysis by finding an equilibrium system that satisfies the 
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boundary conditions and nowhere surmounts the failure 
criterion. The derivation of the failure criterion for wood is 
given here, which is shown by van der Put 1 to be identical 
to a polynomial expansion of the real failure criterion. The 
quadratic polynomial part of this polynomial is shown in the 
following section to represent the critical distortional energy 
criterion of initial yield and thus is not just an expansion, 
but is identical to the real yield criterion. It also is shown 
that the third-degree terms of the polynomial follow the 
theoretical explained mixed-mode I-II equation and repre¬ 
sent special hardening effects by microcrack arrest. The 
fourth-degree and higher-degree polynomial terms have no 
meaning and are zero. As empirical confirmation, the clear 
wood results given by van der Put 1 are discussed together 
with the biaxial data for timber from Hemmer . 2 


Theory 

It was shown for the first time for wood by van der Put 1 that 
the second-degree tensor-polynomial describes initial 
“flow,” which is shown here to represent the extension of 
the critical distortional energy criterion to orthotropic 
materials providing the necessary basis for exact solutions 
according to limit analysis (see Discussion). Because an 
isotropic matrix of a material may sustain large hydrostatic 
pressures without yielding, yield depends on a critical value 
of the distortional energy. This energy W d is found by 
subtracting the energy of the volume change from the 
total strain energy. Thus, for the isotropic matrix material, 
this is 

W d = [ —(ol + ol + x o v + o , + o z o.)+ 

d 1 2E K x - v z> E V x y v 1 z x> 

1+V / 2 2 2 ( 1-2V , yA 

—{r X y + T yz + Tj -1 ——(<L + o y + a z ) \ 

i C 1 ) 
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where a; are the normal matrix stresses; z, the shear stresses; 
E the modulus of elasticity; G the shear modulus; and 
v Poison’s ratio of the matrix material following 2G = El 
(1 + v). 

Wood has to be regarded as a reinforced material and 
initial failure is due to failure of the isotropic matrix. This 
was shown by van der Put, 3 leading to a new fracture 
mechanics theory and a new transformation of the Airy 
stress function. This makes exact solutions possible as 
applied for the derivations of the Wu 4 mixed-mode I-II 
fracture criterion and the derivations of the right fracture 
energies and the relation between mode I and II stress 
intensities and energy release rates. According to van der 
Put, 3 the matrix stresses can be expressed in orthotropic 
stresses as follows. 

The stress in wood a xor is n, times the stress in the matrix 
a x due to the reinforcement in x-direction: a xor = (E x / E)-a x 
= n,-o x , while the reinforcement in v-direction is regarded 
to belong to the matrix, thus c yor = a y and E y = E of the 
matrix. For the shear stress, the multiplying factor is n 6 = 
(2 + v xy + v yx )-G xy / E. Thus, E x , E y , G xy , v xy , and v yx are the 
orthotropic values of wood due to the reinforcements. 

Equation 1 applies for a material with equal tension and 
compression strength. For unequal axial strengths, the 
failure condition, in .v-direction is: (ct v - X)-(a x + X') = 0, 
where X is the tensile strength and -X’ the compression 
strength, as given in Fig. lb. This condition can be written 
as 


X-X'V 

° x 2 J 


X + X' 


or: a x - p x = ±X 


( 2 ) 


(see Fig. 1) and the behavior is identical to that of a material 
with equal tension and compression strengths of X being 
prestressed by stress p x . 

This result follows from the applied linear transforma¬ 
tion. Because Eq. 1 describes a physical property, it should 
be independent of the chosen vector space and according 
to the additivity rule of linear mapping (linear transforma¬ 
tion) is 

f(x + y) = f(x) + f(y), 
or in this case: 


(X-X')/2 



(3) 


f (a - p) = f (tf) + f {-p) giving 
f(a) = f(cT-p)-f(-p) 

Substitution of a x - p x , a y - p r and a, - p z for cr r , o y , and a,, 
respectively in Eq. 1 gives 

f — - — - Oy,or + Py ) + {(?y,or ~Py~ z,or + Pz f + 

V Yi\ yi\ J 


_ „ ^x,or . Px,or . . 

°Z,or-Py -+- I + 


6ZTf ] /n l = 2C(=6EW d /(:\ + v)) 

and after subtraction of f(-p) this is 

2 

6 x,or . 2 . 2 ^x,or^y,or ^ x,or^z,or . 

2 ^ y,or ' ^ z,or ^z,or^y,or ' 


. ^ Px.or | & x or , f 0 , Px.or I _ 

Py+Pz --+ Py~2Pz + - \ a :,or + 


Pz - 2 P y +— k'+ 3 Xrf/n, = ~C p + f(-p) + C 


(4) 


C = 3EW d /(l+v) 


with 


f i~P) = Plor /+ P 2 + Pi - PPx.or /«! “ PPz ~ PzPx.or /«! = C p 


following from inserting a x = p x , a, = p y = /?, and a , = p, in 
Eq. 1. 

Of interest for failure by flat crack propagation is the 
plane stress equation with a, = p, = 0; z xz = z yz = 0; and 
Py = Py,or = P, giving for Eq. 4: 


C'ni 

a 


C'n j 


C’ C'n, 


P-2 


Px,o 


C 


rzi 2p-!-^ L 1 + 3—7 £ tt = 1 


CK 


(5) 


For o y or = z or = 0, Eq. 5 becomes 

x,or . tT x or I r* Px,or | i 

C' n\ C'-n\ P ~ 


This is identical to 


(a x , or -X)(a x , or +X')=0, 
or to 

al or +(X'-X)a x , or -XX' = 0, 
showing that 

(pn 1 -2p xor ) = X'-X, and C'nf= XX' 

The same applies in the perpendicular v-direction for the 
uniaxial tension and compression strengths Y and Y' giving 
C' = YY' and (p xor / n, - 2p) = Y' - Y. This last result is 
to be expected because according to the molecular theory, 
the strength is proportional to the E modulus and thus 
is YY' = XX' / n\ and X' - X = n, (Y - Y'). Then it also is 
P = Px.or / «i = Y - Y' = (X - X')ln,. 

Equation 5 becomes 

2 2 2 
a^__ C> X ' 0 r -O y , or + Oyp L _ O x£! __ ^y z o L+ 3 J ^_ =1 (6) 

C'-nj C' n, C' F C' n, F C' C’-n 2 6 K ’ 
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or 


XX' 

T. 


X 

cr„ 


X' 

a 


- 2 F 12 (T x or C y or + 


YY' 


Y' S 2 


(7) 


where S is the shear strength and 

2 F u = \/C'n ] = l/VXXYF' (8) 

This value of Fu applies for initial stress redistribution and 
microcracking of the matrix and becomes lower by further 
straining reaching a near-zero value at final yield or failure 
initiation. This shows early dissipation of the elastic dis- 
tortional energy, as also is indicated by the orthotropic 
finite element calculation of Gopu. This dissipation of 
distortional energy is according to the incompressibility 
condition and thus follows a minimum energy principle of 
yield. At the end of this stress redistribution, hnal yield 
occurs according to Eq. 7 with F\2 = 0, leading to a com¬ 
bined critical strain energy criterion of the reinforcements 
with a Tresca criterion of the matrix (see Discussion for 
the implications). 


Materials and methods 

The parameter estimation of the first general failure crite¬ 
rion of wood by van der Put 1 was mainly based on uniaxial 
off-axis tension and compression tests on spruce like 
European softwoods and other species from data in the 
literature. 6,7 Because this criterion also applies for uniaxial 
loading, the relations, according to Eq. 18 between the 
parameters, has to be satisfied and it will be shown that 
therefore the uniaxial tests are sufficient for the estimation 
of all parameters. For the empirical verification of this, an 
initial cooperative project was started for biaxial parameter 
estimation, and tests in the longitudinal-tangential plane of 
spruce like European white fir were done by Hemmer 2 on 
torsion tube tests, subjected to tension and internal pres¬ 
sure. However, the verification of Eq. 18 is lacking in litera¬ 
ture and is therefore given here, based on these data. The 
parameter estimation in all cases was based on fitting the 
Tsai-Wu equation, 8 which can be regarded as the tensor- 
polynomial expansion of the real failure equation and has 
the form: 

FjPt + FijOPj + F ijk c i o j o k + .= 1 

(U j, k = 1 , 2 , 3 , 4 , 5 , 6 ) 1 ; 

where a (see Fig. 2) and F are stress and strength tensors. 
In van der Put, 1 it is shown that clear wood can be regarded 
as orthotropic in the main planes and the principal direc¬ 
tions of the strengths are orthogonal (showing the common 
tensor transformations), and third-degree terms represent 
local hardening in stable tests leading for plane stress to 

F 2 g x + F 2 o 2 + F u g 2 + 2F Y2 g 1 g 2 + F 22 g\ + F 66 g 2 6 + 3F 112 g]g 2 + 
3F 221 g 2 2 g 1 + 3F 166 g 1 g 2 6 + 3F 266 g 2 g; + 12F 1266 g 1 g 2 g 2 6 = 1 

( 10 ) 


2 tension is positive 



For less stable tests, failure occurs at yield and Fu is small 
and higher-degree hardening terms are negligible leading 
to Eq. 11, the critical distortional energy equation of initial 
yield derived before: 

F 1 g 1 + F 2 g 2 + F u g 2 + 2F u g 1 g 2 + F 22 g\ + F 66 g 2 6 = 1 (11) 

The relations between the parameters and the requirements 
for parameter estimation by uniaxial off-axis tests is as 
follows. By the uniaxial stress <7„ the stresses in the main 
planes are 

g 1 = g,cos 2 6 G 2 =G,sm 2 6 g 6 = a, cos 0 sin 0 (12) 

Substitution in Eq. 11 gives 

l’\o, cos 2 8 + F 2 g, sin 2 8 + F u g 2 cos 4 6 + 

( 2F 12 + F^g 2 cos 2 0sin 2 6 + F 22 g 2 sin 4 8 = 1 

and expressed in the uniaxial tension and compression 
strengths (X, X', Y, Y') and shear strength S, this becomes 


/ 1 

1 ' 

\ . 2 J1 

i \ 

1 — 


+ G. sm - 8\ — 


lx 

" X 7 , 

I 1 If 

~ Y'J 


2 4 n 2 • 4 

a. cos 6 „ ^ 2 • 2 s\ 2 /-» sm 0 

+-+ 2F 12 g? sin 2 6 ■ cos 2 0 +-+ 

- r, iz r - 


XX’ 
a 2 cos 2 0sin 2 0 


YY' 


(13) 


= 1 


The critical distortional energy equation of initial yield, 
Eq. 13 can be resolved into factors as follows 

f G t cos 2 8 | G, sin 2 8 Af g,cos 2 8 | <7,sin 2 9 | A_ Q 

lx Y j\ X' Y' J ( J 


giving the product of the Hankinson equations for tension 
and for compression, (where X and X' are the strengths 
in the grain direction). This theoretical derivation of the 
empirical Hankinson equations applies according to van der 
Put 1 because 

2 F 12 + 1/S 2 = 1/X'Y + 1/XY’ 


(15) 
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The same can be done with Eq. 10 leading to 

Ffi, cos 2 0 + F 2 o, sin 2 0 + F u af cos 4 0 + ( 2 F n + F 66 ) of cos 2 0 + 
F 22 o 2 sin 4 0 + 3 (F 112 + F 166 ) a, 3 cos 4 0 ■ sin 2 0 + 

3(F 122 + F 266 )<7 ( 3 sin 4 0-cos 2 0 + 12F 1266 ct 4 cos 4 sin 4 0=1 

(16) 


This fourth-degree equation in a, can be written as the 
product of two quadratic equations 


a, cos 2 0 a, sin 2 0 2 . 2 „ 2 n ~ 

+ - L - rz -1 + oysin 0cos 0 C, 




Y 


a, cos 0 a, sin 0 . 2 • 2 „ 2 „ „ , „ 

—-h —-i-1 + cr, sm 0 cos 0-C,, =0 

X' Y' 1 


(17) 



giving the Flankinson equations extended with terms con¬ 
taining hardening parameters C, and C d . For a real failure 
surface the roots will be real and because one of the qua¬ 
dratic equations gives the roots for compression and the 
other for tension, which should be valid for zero values of 
C, and/or C :I as well (reducing then to Eq. 14), the constants 
C, and C d are uniquely determined. 

Performing the multiplication of Eq. 17, the parameters 
of Eq. 16 are known: 

2 F u = l/X'Y + 1/XY’ - 1/S 2 + C,-C d = C t - C d ; 

3(F n2 + F 166 ) = C,/X' + C d /X- (18) 

3 (F 122 + F 266 ) = C,/Y' + C d /Y; and 12F 1266 = C,C rf 

Equation 16 is equal to the empirical criterion of van der 
Put 1 and Flemmer, 2 which was based purely on data fitting 
to the tensor-polynomial, showing that F im and other pos¬ 
sible higher-degree terms are zero or negligible (e.g., F 1266 ). 
This is derived here theoretically, based on the theoretically 
explained hardening terms and Eq. 16 is the exact failure 
criterion of limit analysis. The values C, and C d (depending 
on the stability of the test) follow from fitting the extended 
Hankinson equations for uniaxial off-axis tension and 
compression. 

Theoretically, -^166 — 0 because a 1 is in the plane of the 
crack and is not influenced by the flat crack. For the same 
reason F 122 = 0. Given that F 122 = F l66 = 0, the other constants 
follow from Eq. 18. This will be measurable on perfect 
straight clear specimens. For tension, F 266 = 0 according to 
Eq. 18. This is confirmed by the clear wood measurements 
with n = 2 in Eq. 34 (see Discussion). Timber specimens 
with defects show deviations of the main planes from the 
in-plane stress direction and there are always stresses per¬ 
pendicular to the flat crack plane such that the interaction 
terms F122 and F \66 will not be zero. The high value of 
F \\2 and the similar behavior of F266 indicate that, due to 
the disturbances, random inclined initial cracks outside the 
main planes, are determining showing an extension of the 
in-plane Wu equation to three dimensions. As shown and 
discussed below, it is possible to determine F266 and F\\2 
directly and ^166 and F \22 then follow from Eq. 18. 

As shown in van der Put, 1 F n is very small or zero for 
clear wood. The Flemmer 2 data for biaxial strengths in lon¬ 
gitudinal tension also show for timber that F n is zero at 
yield (see Fig. 3 showing zero slope of the ellipse) and 
thus is 


Fig. 3. Yield criterion for F u = 0 and cr 6 = 0 

1/S 2 = \IX'Y + 1/XY' (19) 

Equation 18 now shows that, due to hardening, F u changes 
from zero at yield to 

2 F u = C,-C d (20) 


Results 

Transverse strength 

In van der Put 1 it is shown that the quadratic polynomial 
like Eq. 11 precisely describes the peculiar flow behavior of 
the transverse compression-tension strength and (rolling) 
shear strength perpendicular to the grain z rol without the 
need for higher-order terms. After some strain hardening, 
the differences between radial-tangential compression 
strength and off-axes compression strength disappear and 
one directionally independent strength value remains. The 
behavior then becomes quasi isotropic and indicates the 
isotropic matrix to be determining. This also applies for 
the second hardening stage after the empty spaces in wood 
are pressed away. For tension perpendicular to the grain, 
only in a rather small region (around 90°, see Fig. 8 of van 
der Put 1 in the radial direction) is the strength higher, 
making it negligible in practice such that a lower bound of 
strength will apply that is independent of the direction. 

Longitudinal strengths 

The longitudinal shear strength in the radial plane increases 
with compression perpendicular to this plane according to 
the coupling term F 266 (direction 2 is radial direction and 
direction 1 is grain direction): 

F 2 o 2 + F 22 g\ + F 66 cj 6 2 + 3F 266< 7 2 ct 6 2 = 1 or 

o* = J (i-q 2 /FHi+cr 2 7n ( 2| ) 

5 V 1 + cg 2 /Y' 

with c = 3F 266 Y'S 2 ~ 0.9 (0.8 to 0.99, see Fig. 4). The value 
of c depends on the stability of the test. When c approaches 
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c ~ 1 (measurements of project A in Fig. 4), Eq. 21 
becomes 



being the Wu equation of fracture mechanics theoretically 
explained by van der Put." 

Equation 21 is thus shown to be an exact equation. As 
derived in van der Put, 3 Eq. 22 does not only apply for 
tension with shear but also for shear with compression a 2 
perpendicular to the flat crack plane. For a high stress <7 2 , 
the crack is closed at o 2 = a c and the crack tip notices only 
the influence of a 2 = o c because for the higher part of cr 2 , 
the load is directly transmitted through the closed crack and 
Eq. 22 becomes 

f= ~ /j(CT f ac) +f^ or a 6 =C + n\a 2 \ (23) 


extension of the Wu equation to the biaxial loading case. 
F\n has the most influence when cr 6 = 0 and when the 
theoretical value F 122 = 0 is reached; then the failure 
equation is 




y 7 


+ ldr + w +3Fin ° 2 ° l=1 (24) 


This can be written as 

o 1 (X'-X) + ot(l + 3F n2 o 2 XX') 

= (1-o 2 /YH1 + o 2 /Y')-XX’ ^ ; 

The critical value of Fm, to just have a closed failure crite¬ 
rion surface, will occur at high absolute values of and <7 2 , 
thus in the neighborhood of (7, == -X'. Inserting this value 
in the smallest term of Eq. 25 gives 

a 2 ( 1 + 3F 112 o 2 XX' + {X'~ X)/(-X')) 

= (l~a 2 /Y)il + a 2 /Y')-XX' 


where a 2 and a c are negative, giving the Coulomb equation 
with an increased shear capacity due to friction: qlcr 2 l. 
However, inserting the measured values of Hemmer, 2 it 
appears that the frictional contribution is very small. The 
microcrack closure stress a c will be about equal to the 
tensile strength: a c ~ -KThe shear strength will be maximal 
raised, at high compression of a c ~ -0.9 F', by a factor 

(l + q(0.9y'-y)/5V2) = (l + 0.3(0.9-5.6-3.7)/9.8-V2) 

= 1.03 

Thus, the combined shear-compression strength is mainly 
determined by an equivalent hardening effect, caused 
by crack arrest in the critical direction by strong layers. 
At higher o 2 stresses, compression plasticity perpendicular 
to the grain (project A of Keenan and Jaeger, 7 see Fig. 
4), or instability of the test (project B of Keenan and 
Jaeger 7 of oblique-grain compression tests) may become 
determining. 

The parabolic strength behavior at biaxial compression 
in the main planes by <7, and ct 2 , with a 6 = 0, determined 
by /•) l2 , given by the longitudinal compression data outside 
the drawn ellipse in Fig. 3, can also be explained by crack 
propagation due to initial off-axis microcracks, giving an 


or 


X' 


{l-o 2 /Y)il + o 2 /Y') 
1 +co 2 /Y' 


where c = 3 F U2 Y'X' 2 


(26) 


Thus, when c = 1, the curve reduces to a parabola and the 
requirement to have a closed failure curve is c < 1. Thus, 

3F U2 <l/Y'X' 2 . (27) 

For longitudinal tension <7 h the right side of Fig. 3 is an 
ellipse, which corresponds to a second-degree equation with 
zero Fll2 and Fl22 due to instability of the test. 

It is not possible to have one precise overall fit for the 
failure criterion for the different failure types. This is 
because the values of the third-degree hardening terms 
depend on the loading state and separate fits are necessary 
for longitudinal tension with zero F\\2 i F l2 , and F\22 and for 
longitudinal compression with hardening, whereby, F\\2 
dominates and Fn is not zero. With the estimates of F 266 and 
F\12 close to their bounds for compression, and with zero 
normal coupling terms for tension, all constants of one 
general failure criterion, Eq. 16, are known according to Eq. 
18, depending on the fitted values of C d and C, from uniaxial 
off-axis tension and compression tests. 



Normal stress in MPa 


Estimation of polynomial constants by uniaxial tests 


Based on data fitting of uniaxial tension and compression 
tests of Hemmer, 2 the values of C d and C, are determined 
here. These are compared with the data and fit of the biaxial 
measurements of Hemmer. 2 

In Fig. 5 of Hemmer, 2 the drawn line and dashed line in 
the figure give the prediction of the uniaxial values based 
on the biaxial measured constants according to the general 
Eq. 16. For comparison, the fits of the extended Hankinson 
equations are also given here. For tension, by scratching the 
nonzero compression term of the product, the extended 
Hankinson equation (Eq. 17) becomes 


cr, cos 2 6 
X 


a, sin 2 6 
+ —- - 


+ <7, 2 sin 2 0cos 2 6 C t = 1 


Fig. 4. Combined shear-tension and shear-compression strengths 


y 


(28) 
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Fig. 5. Uniaxial tension strength and compression strength 


This equation fits the line for tension in Fig. 5 when 
C, = 11.9 IX 2 . The Flankinson power equation (Eq. 34) 
fits in this case for n ~ 1.8 and Eqs. 16, 28, and 34 give 
the same result by only one independent constant n or 
C, 

For compression, the dashed-dot line according to Eq. 
29 is drawn through the measured points at 90° and 15°, 
giving the expected Hankinson power value of n = 2.4 in 
Eq. 34 and C d = 4/(X') 2 in 


cos 2 0 a. sin 2 0 
-+ —- 


V' 


(29) 


The dashed line in Fig. 5, based on the best fit of the biaxial 
measurements of Hemmer, 2 is not a good fit for the oblique 
grain test. This deviation is due to early instability of the 
uniaxial off-grain test. This follows because, for example, 
the ratio of the compression strengths perpendicular to the 
grain and along the grain is 0.2 in the uniaxial tests and 0.1 
in the biaxial tests, showing more hardening in the biaxial 
tests. Thus, strong hardening in the biaxial test will not 
occur under all circumstances and the hardening parame¬ 
ters should be omitted for a consistent safe lower bound 
criterion. 

A fit of the longitudinal tensile data of Hemmer 2 by 
the second-degree polynomial (ellipse) in the principal 
stresses <7, and a 2 (<7 6 = 0) (see Fig. 3) shows F n = 0 and 
F m = F 122 - 0 in the radial plane. Because the Hankinson 
power value for tension n is different from n = 2, there 
must be higher-degree terms for shear (F 166 , F 266 ). The 
strength values according to the data fit in Hemmer 2 are 
(in N/mm 2 ): 

X = 59.5; X' = 41.7; Y = 3.5; Y' = -5.9; S = 9.7 and 
C, = 11.9/X 2 ; C d = 4/ X' 2 

Fl266 can be neglected. It only provides a slight local 
correction. 


The data fit for longitudinal tension with F n = F m = F 122 
= 0, predicts a uniaxial compression strength of X' = 
41.7 MPa, giving the failure at the same hardening state as 
in the oblique grain test (where strong biaxial compression 
hardening does not occur). The constants are 

C, = 11.9 /X 2 = 11.9/59.5 2 = 0.00336; 

C d = 4/X' 2 = 4/41.7 2 = 0.00230 ( ’ 

leading by Eq. 18 to Eq. 32 giving a better fit than the best 
fit of Hemmer, 2 even in the compression range due to the 
application of Eq. 18 and hardening constants. 

In the same way, the data fit for longitudinal compres¬ 
sion, based on the strength values of the best fit of Hemmer, 2 
(for comparison of the results), with c values of c 112 = 0.8 
and c 266 = 0.9 (starting point) leads to Eq. 31. 

Because C„ C d and n values of the Hankinson equations 
are sufficiently close to the published extreme values of n 
for wood with defects, the c values used here also can be 
used in general for all softwoods. For compression, Eq. 16 
becomes 


(Tg 


■[ 1 + 0.9- —+ 0.2- 



(l+^l 

l Y' 

xc 

1 { xj 

l X'J 



1+0.8-^^-0 .77- c^1Cr, 


Y'X' 


X'Y' 


0.41- 


X'Y') 


(31) 


Based on the fit for longitudinal tension (<7, > 0), Eq. 16 
becomes 



(32) 


giving a better overall fit than Eq. 31. Compression 
hardening according to Eq. 31 occurs for low values of a b 
only, and only in the torsion tube test in the radial plane. 
Thus, Eq. 32 more generally represents the failure criterion 
for tension and compression in praxis when n * 2 (in 
Eq. 34). 

For tests and structures, showing early instability and no 
hardening at failure in most situations in practice (causing 
the Hankinson power value to be n = 2), for timber and 
glulam, the determining criterion becomes 



which is identical to Eq. 11 with F u = 0 and IIS 2 = 1 1 X'Y + 
1 IY': 


S 2 X X' XX' Y Y' YY' 


(33) 


Because the virtual work equations, and thus the laws of 
limit analysis, do not apply for hardening, it also is neces¬ 
sary to use Eq. 33 in the codes in all cases for timber and 
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clear wood to replace the now often used, but not valid 
Norris equations. 9 


Discussion 

The failure criterion of wood consists of an orthotropic 
third-degree tensor-polynomial, Eq. 16, first given for 
wood with all aspects by van der Put. 1 This polynomial is 
shown, for the same loading cases, to be identical to the Wu 
mixed-mode I-II equation, 4 which is theoretically explained 
by van der Put 3 in terms of hardening due to hindered 
microcrack extension. Thus, the theoretically explained 
equation, Eq. 16, can be stated to be the real failure crite¬ 
rion for wood. 

The second-degree polynomial part of the failure crite¬ 
rion, Eq. 7 or Eq. 11, is shown to be the orthotropic critical 
distortional energy principle for initial yield showing the 
start of dissipation of distortional energy of the elastic stage, 
and thus is, by this minimal work principle, the exact initial 
yield criterion. After further loading, at final flow, F 12 = 0. 
This means an absence of coupling terms between the 
normal stresses. This is possible only when the reinforce¬ 
ment takes the whole normal loading, causing the matrix to 
fail by shear and the critical distortional energy principle 
thus reduces to the Tresca criterion. The necessary validity 
of the Tresca criterion is confirmed by van der Put 10,11 where 
the strongly increased (sixfold) compression strength under 
the load of locally loaded blocks and the increased embed¬ 
ding strength of dowels is explained by the construction of 
the equivalent slip line field in the specimen using the Tresca 
criterion. In addition, the many apparent contradictions of 
the different investigations are explained by this theory. 
This strong increase of the compression strength is due to 
confined dilatation by real hardening (when the empty 
spaces in wood are pressed away). The Tresca criterion 
satisfies the normality rule and thus inherently the theorems 
of limit analysis for matrix failure. The normality rule thus 
need not apply for the orthotropic criterion. This condition 
is now shown to be replaced by the minimum work condi¬ 
tion for dissipation represented by the yield equation and 
the hardening state constants C d and C, of Eq. 17. 

The distinction between the application of Eq. 10 or Eq. 
11 simply follows from the exponent of the Hankinson 
equation (Eq. 34). Equation 14 shows that the exponent n 
of this Hankinson power formula: 

d, cos" 0 | a, sin" d ^ 

X Y K 

is n = 2 for tension and compression at yield when there are 
no higher-degree polynomial terms. A value of n, different 
from n = 2 thus means the presence of higher-degree terms 
due to hardening after yield given by Eq. 16. 

The initial yield criterion of Eq. 13 or Eq. 14 should 
satisfy both the elastic and the yield conditions at the same 
time. Because the Hankinson power equation with n = 2 
also applies for the off-axis modules of elasticity and because 
this modulus is proportional to the strength, the Hankinson 


power equations with n = 2, Eq. 14, satisfy this requirement. 
Thus, n = 2 is necessary for initial yield. This does not only 
apply for clear wood, but it is mentioned also by Mohler 12 
that n . = 2 for lumber and for glued laminated wood. In 
Kollmann 6 (p. 686, 809), n — 2 is found for clear wood in 
bending and in tension and for the embedment strength on 
(p. 888). The combined shear-compression tests (of Keyl- 
werth in Kollmann, 6 p. 906) show, for off-axis longitudinal 
shear and in the radial plane, that n = 2, showing no harden¬ 
ing terms for the shear strength. This also applies for the 
tangential plane according to the test method of Keenan 
and Jaeger, 7 but in this case not for the radial plane. Early 
instability of uniaxial compression tests and strong harden¬ 
ing in the more stable biaxial compression tests follows 
from the data of Hemmer. 1 

In Eq. 17, the value of the hardening state parameter 
C d for compression or C, for tension, determines 
entirely the hardening state and by that the precise form of 
the measured Hankinson curve. This may also, as an em¬ 
pirical function of the stress level, determine the loading 
curves. 

Because tests in longitudinal compression show other 
hardening rather than tests in tension and thus different 
terms, separate data fits for tension and compression 
are necessary. However, Eq. 32 for tension, also can be 
applied for compression, providing a better general equa¬ 
tion than Eq. 33. Because of possible instability of struc¬ 
tures at yield, hardening can be lacking (as in the oblique 
grain test) and the second-degree polynomial, or critical 
strain energy yield condition (Eq. 33), applies as the ulti¬ 
mate failure condition for the codes. 


Conclusions 

A derivation is given of the orthotropic extension of the 
critical distortional energy principle of yield, showing energy 
dissipation according to a minimum energy principle. This 
leads at final yield to an extended Tresca criterion for matrix 
failure. For stable tests, hardening is possible, leading to a 
general polynomial failure criterion of the fourth degree. 
This criterion can be resolved in two factors: for com¬ 
pression and for tension leading to extended Hankinson 
equations for uniaxial loading, making determination of all 
constants of the general failure criterion possible from 
simple uniaxial oblique-grain compression and tension tests 
instead of problematic biaxial tests. The extended Hankin¬ 
son equations also show necessary relations between the 
constants of the general criterion, providing necessary data- 
fitting conditions. For compression and for tension, the 
value of one hardening state constant fully describes the 
form of the hardening curve at different hardening stages, 
providing the possibility of a new hardening description. 
The normality rule does not apply for orthotropic wood 
because the matrix is determining initially. Because dissipa¬ 
tion at yield is according to a minimum energy principle, 
this rule can be replaced by the hypothesis of intrinsic 
minimum work. 
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